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ABSTRACT 


In spite of the fact that it was predicted more than four 
decades ago, the axion is still the most convincing answer to 
the strong-CP problem and a well-motivated dark matter can- 
didate, sparking a variety of elegant and ultra sensitive exper- 
iments based on axion—photon mixing. This article examines 
the first occurrences of non-one-dimensional scattering and de- 
rives enhanced probability for a specific model. We selected 
to characterise the inhomogenous magnetic field surrounding 
the solenoids axis by fitting a Gaussian distribution around its 
centre and calculating the three-dimensional axion conversion 
cross section from photon scattering for this particular sce- 
nario. An analysis and simulation in Python were performed 
on the correlation between the transmitted axion cross sec- 
tion and arbitrary variables.Furthermore, we talk about the 
different approaches that have been suggested to look for dark 
matter and solar axions. As a result of their application and 
development, the concept of an invisible axion has come un- 
der considerable pressure. These findings compel us to take 
a closer look at how experiments intended to discover axions 
are carried out, and to understand how to use this research 
on a large scale to understand the concepts of Dark matter 
detection. 


KEY WORDS: Axion; primakoff effect; Photon—axion conversion; ; strong- 
CP problem; Dark matter; 
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1.1 STRONG CP PROBLEM 


The Standard Model (SM) of elementary particles rose to prominence in the 1970s as 
the most accurate representation of all fundamental interactions other than gravity. 
It has since proven to be extremely successful, explaining virtually all important data 
with only a few parameters. However, it was recognised early on as a puzzle: one would 
not anticipate the strong interactions within the SM to retain parity P or the product 


CP of charge conjugation C with parity. 


P and CP are observed to be conserved by strong interactions and electromagnetic 
interactions. On the other hand, weak interactions violate P, C, and CP. The problem 
with the SM is that unless an unforeseen cancellation occurs, the P and CP violation 
of the weak interactions causes P and CP violation in the strong interactions. The 


Strong CP Problem is what this is known as. 


1.1. STRONG CP PROBLEM 


1.1.1 @acp and Neutron dipole moment 


A parameter, Oacp, which appears as the coefficient of a P and CP odd term in the 


action density, controls the amount of P and CP violation in strong interactions. 
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g ny. 
Lym =... +0 $§_G? Gov sie 
SM + QCD 359 Tw (1.1) 





where the G¥,a = 1,2,...8, are the Quantum Chromodynamics (QCD) field strengths 
(QCD), Gow = sehr Ge ,, and g, is the QCD coupling constant. We use units with 
h = c = 1 and conventions with the Minkowski metric (1) = diag(+1, —1, —1, -1) 


0123 — +1 unless otherwise stated. All other terms in the SM action density, 


and € 
i.e. the terms that contribute to its various achievements, are shown by the dots. 
Eq. (1.1) shows the one term that is not a success. O@gcp is an angle, i.e. it is 
cyclic with period 27. QCD depends on #gcp because of the existence in that theory 
of quantum tunneling events, called ”instantons”, which violate P and CP if 6acp 
differs from zero or 7. Since in actuality the strong interactions obey P and CP, as 
well as can be observed, Oacp must be approaching one of the CP conserving values. 
The experimental upper limit on the neutron electric dipole moment provides the best 
constraint: |d,| < 3-107°° e cm( 90% CL). For small @gcp the contribution of the 
term shown in Eq. (1.1) to the neutron electric dipole moment is of order 


e 1 
ane 7 ~3- iC as bacp ecm (1.2) 





dn ~ Pgep My + ma AQcd Mn 
where m, and mg are the up and down quark masses, m, is the neutron mass, and 
Agcp is the QCD scale. 9gcp should therefore be less than of order 107° (modrz) . 
$acp = 0 or 7 is unexpected in the SM because P and CP are violated by the weak 
interactions. CP violation is introduced by giving apparently random phases to the 
Yukawa couplings that give rise to the quark masses. The total phase of the quark 
mass matrix flows into Oacp, which is of order one in general. The Strong CP Problem 


is the puzzle of why 6gcp, which should have been of order one, is less than 107°. 


1.1. STRONG CP PROBLEM 


1.1.2 The PQ solution 


Peccei and Quinn (PQ) proposed a modification of the SM that provides a solution 
soon after the Strong CP Problem was identified. They postulated a Upg(1) symmetry 
that 


1) has the classical action’s exact symmetry, 
2) is spontaneously broken, and 


3) has a colour anomaly 


i.e. it is explicitly broken by the non-perturbative QCD instanton effects that make 
physics depend on the value of @gcp.When this recipe is followed, the axion decay 
constant is substituted by (x)/f¢, o(x) with x = (x°, x’, 2?, x*) is a dynamical pseudo- 
scalar field and fy is a quantity with is an energy-dimensioned quantity. The vacuum 
expectation value of order f spontaneously breaches Upg(1) symmetry. The associated 
Nambu-Goldstone boson is ¢(xz). The non-perturbative instanton effects that make 
physics rely on @gcp introduce an effective potential for ¢(x), according to Weinberg 
and Wilczek (WW). The minimum of this effective potential was later shown to be at 
(x). The Strong CP problem is solved after the (x) field settles there. 


The PQ mechanism alters the SM’s low-energy effective theory by introducing a light 
pseudo-scalar particle known as the ”axion,” which is the quantum of the ¢(2) field. 
The Axion’s properties are largely determined by the value of the axion decay constant 
fs. The axion mass mg is inversely proportional to fz, as are all of its interaction 
strengths. fy is of order the electroweak scale in the original PQWW model, implying 
an axion that is relatively strongly coupled and heavy, i.e. mg of order 100 keV. Several 
laboratory experiments, including unsuccessful searches for axions in beam dumps and 
rare particle decays such as Kt — a* +a, as well as stellar evolution constraints, 
quickly ruled out the PQWW model. The latter is due to the fact that stars emit 
weakly coupled axions from their cores but only photons from their surfaces. If axions 
exist, stars have an additional mechanism for losing energy, causing them to evolve more 
quickly. Axion models with fs < 10° GeV are ruled out when the negative results from 


accelerator-based axion searches are combined with stellar evolution constraints. 


The original PQWW model is unworkable, yet the underlying idea of Peccei-Quinn 


symmetry and its accompanying axion is still valid today. Upg(1) does not need to be 


1.1. STRONG CP PROBLEM 


broken at the electroweak scale, according to Jihn E. Kim and colleagues. It could be 
broken at any energy level, such as the hypothetical ” grand unification scale” of 1015 
GeV. When fy is large enough, the axion is exceedingly light and very weakly coupled 
(mg ~6-10~° eV for fy = 10" GeV)) : all axion production and interaction rates are 
suppressed by about 25 orders of magnitude compared to the PQWW axion. The 
“Invisible axion” was created by arbitrarily increasing fy, a solution to the Strong CP 
Problem that elegantly avoids the limits imposed by laboratory searches and evolution 


of star. 


Thankfully, cosmology saved the day. Indeed, axion field oscillations must begin suffi- 
ciently early in the history of the cosmos for ¢(x) to relax to zero. Since the oscillation 
period is 2/mg, the axion must be suitably heavy. Because of the universe’s finite age, 


there is a limit to how small mg, or equivalently how large fi, can be. 


CHAPTER 2 


AXION - A DARK MATTER 
CANDIDATE 
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2.1 AXION and DARK MATTER 


Several cosmological parameters have been measured with high accuracy recently, and 
the results have shown an universe of which only a small portion can be viewed di- 
rectly. Deuterium abundances anticipated by the hypothesis of big bang nucleosynthe- 
sis (BBN) have been measured, and it has been shown that baryonic matter accounts 
for just 4% of the total energy density of the universe . According to evidence gathered 
from the cosmic microwave background, along with supernovae searches, galaxy surveys 
and other measurements, the great bulk of the universe is composed of gravitational 
”Dark Matter” (26 %) and negative pressure ” Dark Energy” (70 %). 


The presence of ” Dark Energy” has just recently been demonstrated (mainly by cos- 
mic supernovae surveys conducted in the last decade), although the existence of ” Dark 
Matter” has been known since the early 1930s. That is when Fritz Zwicky observed 


that the member galaxies of the Coma cluster were moving much too fast to be gravi- 


2.1. AXION and DARK MATTER 


tationally bound by the bright matter, which led to the discovery of the Coma cluster 
itself. The two possibilities are that they were unbound, which would indicate that the 
cluster should have split apart billions of years ago, or that there was a significant quan- 
tity of undetected ”dark matter” holding the system together. Because of these early 
discoveries, evidence for dark matter has been accumulating on sizes ranging from the 


size of dwarf galaxies (kiloparsecs) to that of the whole visible universe (gigaparsecs). 


At this time, the best dark matter possibilities appear to be unknown non-baryonic 
particles that were left over after the big bang, according to current research. In terms 
of interactions with standard-model particles such as baryons, leptons, and photons, 
they would be considered to have only the weakest interactions. In light of recent 
research into the creation of cosmic structures, it appears that the vast majority of dark 
matter is ”cold,” that is, non-relativistic, at the time of galaxy formation’s inception. 
Relativistic dark matter would likely to stream out of early density disturbances since 
it is collisionless, essentially smoothing out the cosmos before galaxies had a chance to 
form. The galaxies that we view now are typically surrounded by enormous halos of 
dark matter that stretch far beyond than their luminous limits, a phenomenon known 
as haloing. Dark matter density near the solar system is estimated to be about 09 pM 
~ 0.45 GeV cm~? based on measurements of the Milky Way’s rotation curves (together 


with other observables such as microlensing surveys). 


There are two main types of dark matter candidates: the general class of Weakly 
Interacting Massive Particles (WIMPs), which includes the supersymmetric neutralino 
as an example, and the axion, which was predicted as a solution to the so-called 
”Strong CP” problem but has not yet been seen. Despite the fact that both particles 
have strong motivations, this explanation will concentrate solely on the axion. Axion 
is a light chargeless pseudo-scalar boson (negative parity, spin-zero particle) predicted 
by the violation of the PQ symmetry, which is explained in detail in Chapter 1. In the 
late 1970s, this symmetry was first proposed as a way to explain why charge (C) and 
parity (P) appear to be preserved in strong interactions, despite the fact that the QCD 
Lagrangian contains an overtly CP-violating component. Although this CP-violating 
term should have resulted in a readily observable electric dipole moment in the neutron, 
no such moment has been found to extremely high precision in any of the neutron’s 


experiments. 


The spontaneous symmetry breaking (SSB) scale of the Peccei-Quinn symmetry fy, 
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is the essential parameter dictating most of the axion’s features. The mass and axion 


coupling are both inversely proportional to fy., with the mass specified as 


10'? GeV 
ip POR RY (a=) (2.1) 
fe 
and the axion photon coupling constant (g¢,) expressed as 
a 


where C’ is a dimensionless model dependent coupling parameter and ais the fine 
structure constant. In general, C' is assumed to be ~ 0.97 for the KSVZ (Kim- 
Shifman-Vainshtein-Zakharov) class of axions and ~ —0.36 for the more pessimistic 
grand-unification-theory inspired DFSZ (Dine-Fischler-Srednicki-Zhitnitshii) models. 
Because interactions are proportional to the square of the couplings, these C’ values 
tend to limit the feasible axion-to-photon conversion rates to a few orders of magnitude 


at any given mass. 


Initially, f was thought to be on the electroweak scale (fy ~ 250 GeV), implying an 


axion mass of order 100 KeV and strong couplings visible in accelerators . 


Searches for axions in particle and nuclear experiments, together with astrophysics 
limitations, quickly dropped its possible mass to my < 3 x 107? eV, equating to fy > 
10° GeV. 


These low-mass axions were first assumed to be undetectable since their couplings are 


inversely proportional to fy. They were dubbed ” invisible” axions. 


It was discovered that a general lower limit on the axion mass may also be derived 
from cosmology. Axions would be abundant at the moment of the big bang, thanks to 
several mechanisms which is discussed in this paper. The overall contributions to the 


universe’s energy density from axions formed using the vacuum misalignment approach 


7/6 
Oy (=) (2.3) 


Me 


may thus be written as 


This implies that the axion mass is mg > 10~° eV (any lighter and the axions would 
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overclose the universe, 4 > 1). When the astrophysical and experimental constraints 
are added together, the axion has a mass range of weV — meV, with the lower values 
being more likely if the axion is the primary component of dark matter. The axions 
produced in the early cosmos around the QCD phase transition, when the axion mass 
switches on, would have momenta. The temperature of the surrounding plasma was 
T ~ 1 GeV, while the temperature of the surrounding plasma was ~ 107° eVc7!. 
Furthermore, because such axions are so weakly coupled, they will never be in thermal 
equilibrium with anything else. This means they’d be non-relativistic ”cold” dark 
matter right away, and they’d be able to build structures around density perturbations 


rather soon. 


A massive halo of particles travelling at relative velocities of the order of 10~* c would 
make up the axion dark matter in the galaxy today. It’s unknown if any or all of the 
axions would be gravitationally thermalized, but they’d have to to be going slower than 
the local escape velocity of 2 x 10-* c to be bound in the galaxy. Non-thermalized 
axions may still be oscillating into and out of the galaxy’s gravitational well. These 
axions would have extremely small velocity dispersions (on the order of 107!” c), and 
the variances in velocity from different infalls (first time falling into the galaxy, first 
time flying out, second time falling in, etc.) would be correlated with the galaxy’s 


development. 


2.2 AXION MODELS 


Axion models known as ‘Ivisible Axion Models’ because of their extremely weak cou- 
plings are still possible. The PQ symmetry is decoupled from the electroweak scale in 
an invisible axion model, and it spontaneously breaks at a considerably higher tem- 
perature, lowering the axion mass and coupling strength.The KSVZ and DFSZ models 
are two benchmark, invisible axion models. Both of these theories produce an axion 


with permissible mass and couplings. 


The standard Higgs doublet is the only Higgs doublet in the KSVZ model. An ad- 
ditional electroweak singlet scalar field introduces the axion as its phase. The known 
quarks cannot couple directly to such a field because this would result in unacceptably 


large quark masses. Instead, the scalar is connected to another heavy quark, this time 


2.2. AXION MODELS 


an electroweak singlet. The interactions of the heavy quark with the other fields then 


induce the axion couplings. 


The DFSZ model, like the PQWW model, has two Higgs doublets and an additional 
electroweak singlet scalar. At the PQ symmetry-breaking scale, it is the electroweak 
singlet that acquires a vacuum expectation value. The scalar couples to quarks and 


leptons through interactions with the two Higgs doublets, rather than directly. 


String compactifications produce PQ symmetries, which are always broken by an in- 
stanton. While this might seem to make the axion an excellent dark matter candidate, 
string models favour a PQ scale considerably higher than what cosmology allows, pro- 
vides an assessment on the current scenario. It is difficult to lower the PQ scale below 
1.1 x 10'© GeV, as mentioned in ’Axions in string theory’ by Svrcek and Witten , while 


it is much easy to raise its value. 


2.2.1 Electromagnetic coupling 


Axion couplings to other particles are inversely proportional to f, in general, but the 
exact strength of these couplings varies depending on the model. The coupling between 


axions and photons, for example, can be represented as 
Loy = Jory 9E-B (2.4) 


E and B are the components of the electromagnetic field. The coupling constant is 


Gy 
Gory = (2.5) 
ory whe 


where g, is a constant containing the model dependence, 10° GeV < fg < 10° GeV is 
the axion decay constant, of the order of the PQ scale. Explicitly, 


m3 (H- aes), 20 


where z is the mass ratio of the up and down quarks , N is the axion colour anomaly, 








and E is the axion electromagnetic anomaly. The phrase holding the mass ratios of 


light quarks has a value of 1.95. E/N is the ratio that causes model dependence. E and 
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N are connected in grandunifiable models, and E/N = 8/3. This is where the DFSZ 
axion model comes in, with g, = 0.36 in this example. E = 0 and g, = 0.97 for a 
KSVZ axion. 


The development of the KSVZ and DFSZ axion models demonstrates that an axion 
can address the strong CP problem. While significant in terms of itself, the axion also 


offers an intriguing possibility for the universe’s cold dark matter. 


2.3 COSMOLOGICAL AXION PRODUCTION 


2.3.1 Axion dark matter properties 


Axions meet both of the requirements for cold dark matter: (i) a non-relativistic pop- 
ulation of axions could exist in sufficient quantities in our universe to provide the 
required dark matter energy density, and (ii) they are effectively collisionless, with 
gravitational interactions being the only significant long-range interactions. 


(2.7) 


102 GeV 
mg ~ 6 x 10-SeV oS) 


fe 


Axion dark matter is non-relativistic, despite its small mass, since cold populations 
are formed out of balance. Cold axions are created by three mechanisms: vacuum 
realignment, string decay, and domain wall decay. The history of the axion field as the 
cosmos expands and cools is discussed in this section to see how and when axions are 
formed. We also go over vacuum realignment production in brief, note here that this 
mechanism will always contribute to cold axion populations, and it may be the sole 


one, as mentioned below. 


2.3.2 Topological production of axions 


There are two crucial scales in the generation of dark matter axions. The first is the 
temperature at which the PQ symmetry breaks, denoted by the symbol Tpg. Which 
of the three methods contributes the most to the cold axion population is determined 


by whether this temperature is higher or lower than the inflationary reheating tem- 


10 
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perature, Tp. The second scale is the temperature at which the axion mass becomes 
important due to non-perturbative QCD effects. The QCD effects are insignificant 
at high temperatures, and the axion mass is negligible. At a critical point, t;, when 
met; ~ 1, the axion mass becomes essential. 7, ~ 1 GeV is the temperature of the 


cosmos at ¢1. 


At early dates and temperatures larger than Tpg, the PQ symmetry remains unbroken. 
It spontaneously breaks at T’pg, and the axion field, proportionate to the phase of the 
complex scalar field obtaining a vacuum expectation value, can have any value. The 
phase changes continuously from one horizon to the next, altering by order one. Axion 


strings emerge as topological flaws. 


If Tpq > Tr, the axion field is homogenised over enormous distances and the string 
density is diminished by inflation to the point where our visible universe is exceedingly 
unlikely to contain any axion strings. The axion field is not homogenised in the scenario 
Tpq < Tp, and strings radiate cold, massless axions until non-perturbative QCD effects 


become substantial at temperature, 7}. 


There is no agreement on the expected spectrum of axions from string radiation, and 
there are two options. Strings either fluctuate several times before entirely decaying 
and axion production is sharply peaked around a dominant mode, or much faster 
decay occurs, generating a spectrum inversely proportional to momentum. Rapid decay 
generates approximately 70 times less axions than slow string decay, resulting in distinct 


cosmic constraints on axion mass. 


The axion strings become the limits of N domain walls as the cosmos cools to T;. When 
N = 1, the walls quickly emit cold axions and disintegrate (domain wall decay). Because 
the vacuum is multiply degenerate and there is at least one domain wall per horizon 
if N > 1, the domain wall problem occurs. These walls will eventually dominate the 
energy density, causing the cosmos to expand as S « t?, where S denotes the scale 
factor. Although different solutions to the domain wall problem have been presented 
by S. Chang, C. Hagmann, and P. Sikivie (Phys. Rev., D59:023505, 1999). we will 


assume N = 1 or Tpg > Tp in this case. 


As a result, if Tpa < Tr, string and wall degradation contribute to axion energy density. 
If Tpg < Tp and the axion string density is diluted by inflation, these methods have no 


meaningful contribution to the density of cold axions. Then, only vacuum realignment 
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will make a major contribution. 


2.3.3. Vaccum realignment method 


Cold axions will be created independently of Tz via vacuum realignment. The details of 
this method are out of scope of this paper, but the basic mechanism is as follows. The 
axion field amplitude at Tpg can be any magnitude. If Tpg > Tp, inflation will cause 
homogeneity, and the axion field will be single valued over our observable universe. 
Non-perturbative QCD effects generate an axion field potential. When these effects 
become significant, the axion field’s potential begins to oscillate. As non relativistic 
matter, these oscillations do not decay and add to the local energy density. Thus, 
independent of the inflationary reheating temperature,vacuum realignment produces a 
cold axion population. 

( for further study, please see the article: arXiv:1904.05707 |astro-ph.CO]) 
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3.1 INTRODUCTION 


Low-mass axions were once assumed to be ” invisible” to typical observational technolo- 
gies because to their extremely slow decay rates and weak interactions with hadronic 
matter and electromagnetism. Sikivie, on the other hand, demonstrated that the in- 
verse Primakoff effect substantially accelerates the decomposition of dark matter axions 


under a strong static magnetic field. 


One photon is ”replaced” by a virtual photon in a static magnetic field, while the 
other maintains the axion’s energy, which is equal to the rest-mass energy (mgc”) plus 
the nonrelativistic kinetic energy. In Eq. (2.4), B is essentially converted to Bo, the 
static magnetic field. As B?, the axion’s decay rate effectively increases. The Feynman 


diagrams of the axion photon interaction for the two scenarios are shown in Figure 3.1. 


Sikivie presented a method for detecting axion dark matter based on the Primakoff 

effect, in which a microwave cavity permeated by a strong magnetic field is utilised 

to resonantly increase the number of photons produced by axion decay. When the 
MgC? 


resonant frequency f * ——, where h is Planck’s constant, the axion-photon conversion 
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Y Y 


a-->- a-->- 


B, 
(a) (b) 


Figure 3.1 Axion decay into photons as illustrated by Feynman diagram . a) The 
vacuum conversion. b) In a static magnetic field, the inverse Primakoff effect 


is enhanced. There is also a slight correction owing to the axion’s kinetic energy, but 
it is negligible for cold dark matter (42 ~) Liye he This axion dark matter detector is 


known as a haloscope, according to Sikivie. 


The two-photon coupling of axions or axion-like particles allows for transitions be- 
tween them and photons in external electric or magnetic fields as shown in Fig.3.1. 
This mechanism serves as the basis for the experimental searches for galactic dark 
matter axions and solar axions . The astrophysical implications of this mechanism 
have also been widely investigated and reviewed . The phenomenological consequences 
of an extremely light or massless axion would be particularly interesting in several as- 
trophysical circumstances such as polarization of radio-galaxies and quasars, the diffuse 


X-ray background, or ultra-high energy cosmic rays. 


One fascinating cosmological consequence of this technique is photonaxion conversion 
induced by intergalactic magnetic fields, which results in the dimming of distant light 
sources, particularly type Ia supernovae, which are employed as cosmic distance indi- 
cators. SNe Ia at redshifts of 0.8 S z S 1.7 appear fainter than expected from the 
luminosity-redshift relation in a decelerating Universe,a finding usually interpreted as 
evidence for acceleration of the cosmic expansion rate and thus for a cosmic equation 


of state (EoS) dominated by a cosmological constant. 


A slowly changing scalar field, or something more stranger. The darkening generated 
by photon-axion conversion could imitate this phenomenon, offering an alternate ex- 
planation to cosmic acceleration. Although it still required some non-standard fluid 


to account for the Universe’s flatness, this model appeared capable of explaining the 
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SN-dimming via a fundamentally different mechanism. 


However, if light from distant SNe Ia is partially transformed to axion-like particles and 
reaches us, the same mechanism would apply to all remote sources of electromagnetic 
radiation. As aresult, it appears important to update the many arguments constraining 
photon-axion conversion in intergalactic magnetic fields, particularly the limitations 
originating from cosmic microwave background (CMB) spectral distortions and quasar 


(QSO) spectral dispersion. 


3.2 LAGRANGIAN AND EQUATIONS OF MO- 
TION 


This chapter examines the conversion of axions to photons in the lack of cavities or 
reflecting walls for the photons in a static magnetic field. 
In this paper we will consider a system described by the action for the axion and 


electromagnetic fields:, 


1 1 ~ 
oS i da[-p FY Fy — 5 (OudO"o + mo?) + SOF” Fi] (3:1) 


with 


J 1 Va 
Fw = 5" P Fag (32) 


where @ is a real scalar field representing the axion, F),, is the electromagnetic field 
tensor, and g is the coupling constant. 

We are interested in the physical condition of photons and axions propagating in the 
presence of a static magnetic field as a background. We will assume that this magnetic 
background B is dependent on r and points in any direction denoted by n, such that 


the important part of the interaction term is 


i: d‘xBoEn (3.3) 
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where 6 = gB(r) and Ex is the component of the photon electric field parallel to the 
background magnetic field 
Using the Coulomb gauge and ignoring photon field components that do not couple to 


the axion via the magnetic field, (3.1) becomes 


S= 5 | d‘z[0, AO“ A + (0,606 + mo”) — 2860,A] (3.4) 


where A is a magnetic vector potential. Choosing the temporal gauge Ag = 0, then 


the relevant electric field component becomes Ey = 0,A 


(0,0" + m*)o(F, t) = B(F)OA(F, t) (3.5) 
(Q,O")A(F, t) = —B(r) 20, t) (3.6) 
where = (0,1,2,3) Now considering the scenario of an incident wave with time 


dependency f (re, where f (7) represents the fields ¢(7,t) or A(7,t). In free space 


regions where the magnetic field vanishes (zones I, HI of figure 1) the f(r) for A is 
bike? 








etike-® and for @ is ee" where kp is photon wave vector, kp is axion wave vector. 


We write the fields as follows after performing a perturbation expansion in powers of 
£. 
(Ft) = 6% (Ft) +46 (Ft) +... (3.7) 


A(7, t) = AO, t) + AD (Ft) +... (3.8) 


where the equations of motion (3.5) and (3.6) becomes 


(0,,0" +m?) oO (F, t) = 0 (3.9) 
(0,0") AM (7, t) = 0 (3.10) 
(0,0" +:m?)o (F, t) = B(7)0,AM (F, t) (3.11) 


3.3. GREENS FUNCTION CALCULATION 


(G,")AM (Ft) = —B(F)H,0 (Ft) (3.12) 


In the equations (3.11) and (3.12), we can see that the solutions of the unperturbed 
fields (3.9) and (3.10) becomes the sources for the perturbed fields. 


3.3 GREENS FUNCTION CALCULATION 


Calculating the solution for the equation (3.11) using Greens’s function method, and 


we are interested in the solutions for G,,,(7, t) 


(0,0" + m2)Gn(F,t) = —53(F)4(t) (3.13) 


A direct calculation in fourier space gives us the desired Green’s function, 
where Gin(r,t) = 





ae 2 Be gD 2 
poe Writing k* — m* =w* — wy 





1 i(k.F—wt) 
Gn(Ft) = / Pk; (3.14) 
with w, = Vk? +m?. 


In order to deal with the singularities in carrying out the w integration, we will treat w 
as a complex number and treat the integration as a contour integral running along the 
real w axis. The singularities on the contour are avoided by pushing the singularities 
at w = Ww, infinitesimally away from the real axis. Since we want the retarded Green’s 
function solution, that is G,,(r,t) = 0 for t = 0. we integrate on the real w axis, with 
poles +i; shifted to +w, —7te. When t < 0, closing the contour up gives 0; when t > 0, 


closing the contour down gives the retarded Green function 


> 


Get) = Torres: / = el F gin (wt) (3.15) 


For our case of incident photons from the left, we take the incident photon field as a 





monochromatic wave A = ei(*»-?-t) and 6 =0, 
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so that the inhomogeneous equation (3.11) reads 

(0,0" + m)o (#, t) = —i@B(F) ete? #9 (3.16) 
then particular solution construction using the Green’s function for the equation (3.11) 


will be given by 


=. 


oF, t) = —i@ i rid’ Gr —#',t—t)B(7') be") (3.17) 


=-ia f dr'a@ yet [ avak@-F't-t)e™ — (a.18) 
V = 


CO 
where we have assumed (7°) = 0 outside the interval 0 to L. Substituting the Retarded 


green’s function (3.15) in the above equation we get 





Ase 2 ds oo 1 
a 3,,/ DI\ Kikp.? aa _ 4 
= io f abr Br \e / dt’ — O(t Ome 





(3.19) 
~ ak (ot ol oa) 
x i wey ek (FF )sin ww, (t =) Hew 
—co Wk 
— Toor 2 —1 d’k Riso 
= ine f Br! BF") elke? | ; elk FF ) 
V co (277)? wy, (3.20) 


x / dt! o(t — t' et) sin w,(t — t') 


the integral over t’ can be done by subtituting T = t—t’ , there are only contrubutions 
from T > 0 because of the 6(T). after doing the integral on t’ the equation (3.20) 


becomes, 


a fo ag. Py eS ak Too ol Ww) 
= 22 yt 307 SI\ otkp.F ik.(7—7') fe S, 
= —1IWwe [a r B(r Je ia (Qn)3 = e x we = aD (3:21) 


writing |(r — 7 ’)| = R and now solving for the volume integral over k , the above 


equation becomes 


- Eee, ol © kdk eof 
= ine f dr! B(F’) elk? x i i elk cost tog sd 
y Ore, wa I 
(3.22) 
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— 7. tl 2 1 =, } 
= ine i Br! BF") x i, BaintkE) (3.03) 
V 0 





(2n)2R we — w? 
solving separately for the integral over k after substituting w2 = k? + m? 
~ k sin (kR) 
5 5S dk (3.24) 
9 k+m*-—w 


The integral on k can be performed as usual in the complex plane by closing the contour 
in Upper Half plane for r—r’ > 0 and for r—r’ < 0 we close the contour in the Lower 


Half Plane , which gives the results, 


oS ES Lee ihe pan SO 
0 (k+ Va? —m2)(k — Vw? — m2) 2 e-WVP MPR. ppl <0 


(3.25) 
substituting the result (3.25) in to equation (3.23) thus we get the perturbed first order 


transmitted and reflected wave function. 


Prrananited Ft) = ge i dr! O(r — r')B(F yer 
Vv 





5 of eniVe me? (FF) 
6) ected, t) = —e ™ / Br! O(r! — r)B(F ei *"| £ 
V 


3.3.1 Asymptotic limit of the wave function 


We will now demonstrate that for high values of r, the expression (3.26) will now be 
reduced. For ease of calculation we will substitute « = /w?—m?. In a scattering 
experiment, because the detector is located at distances (away from the target) that 


are significantly greater than the size of the target (Figure 3.1), we have r > 1’, 
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Detector 







Target area 


>y 


Figure 3.1. The distance r between the target and the detector is too large in 
comparison to the target’s size r’. ie: r >> 1’ 
Image courtesy : Zettili, Nouredine. Quantum Mechanics: Concepts and Applications. 
Chichester, U.K: Wiley, 2009. Print. 


where r represents the distance from the target to the detector and r’ represents the 


size of the detector (see Figure 3.1). If r >> r’, we can approximate « |r’—7’| and 
|F—F"|~ by 


P| = bP? = OFFER? & wr — oP = er RF" (3.28) 
r 





t ft 7 Lf, FY 1 
ee ne 3.29 
FF) ril—F- PF? > ( ar ) : ee) 


substituting (3.28) and (3.29) into the equations (3.26) and (3.27) and simplifying, we 


arrive at the equations for Oat t) and Go mmitea(® i), 





i(wr—t) [ 
1 e€ Ww 
G snited Ce t) = 


re [ero — r')B(F Vet aa (3.30) 


r 


Oo ja — 2 [ ero! —na@ ye! ad (3.31) 
V 
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CHAPTER 4 


CROSS SECTION CALCULATION 
WITH DIFFERENT MAGNETIC FIELD 


Contents 











4.1 FINITE SIZED SOLENOIDAL GENERATED 
POTENTIALS 


4.1.1 Gaussian Distributed Magnetic Field 


We want to be able to acquire quantifiable values that can be used in a laboratory 
experiment, therefore we must explore a more realistic function to describe the magnetic 
field created by the solenoid in order to achieve this. As an initial model, we chose 
to characterise the inhomogeneous magnetic field around the solenoid’s axis using a 


Gaussian distribution around its centre. 


2 


=f 


B(r) = Boe 2 (4.1) 





From equation (3.30) , the scattering amplitude f(0,@) of the transmitted wave func- 


tion Go mnitea "> t) is given by 


f(0,¢) =— | Br’ A(r' — r) BP el Rot 8) (4.2) 
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) = g Bir '). Where 


where for simplification we substitute ¢ = k, + & and 6( 
= r! 2 
B(r') = Boe ® , then our above equation becomes 
iW 3) yl 2 = 
100; = le ies O(r—r')gBoe ® ef (4.3) 
T 
Evaluating the volume integral using cylindrical coordinates 
a ee cae = ee te 
— By | cerar' dge'"” saa dz (4.4) 
ar 0 0 0 
(4.5) 


fore) /2 20 
= 980 | cera! | doe’ res? 
0 0 


=A 
Using the integral definition for Bessel function of First Kind 
Jo(qr’) = Nie det” © we have 
wo Sct 
=> 9B | ew r'dr’ Jo(qr’) (4.6) 
0 
wl R qe R? 
= gBp — Xx —e 7 4.7 
9gP0 5} 5} (4.7) 
(4.8) 


Therefore our scattering amplitude 
GT R* is 


f(@,¢) = ae gBoR’ e 


A: 
(4.9) 


where the explicit dependence of q on the scattering angle @ is given by 


q’ = 2K?(1 — cos 0) = 4K? sin*(6/2) 


(4.10) 


where K = k, which is the incident photon momentum. The differential cross section 
is related to the scattering amplitude 


do 
= (6,6) = |F(6,0)0 
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The probability density for detecting the scattered projectile at a given angle is the 
straightforward interpretation. 


Therefore the differential cross section is 


= 2 
do wl 2p? 
—(6,¢) = | —gBoR? | x e (4.11) 
dQ 
Integrating the differential cross section oe across the complete solid angle ( 47 stera- 


dians) yields the total cross section o: 


20 
Ctot = a0 = / [3 — sin ddédy (4.12) 
an 


Finding Total cross section from equation (4.11) 


27 @L 5 2 ee? 
C= we BR e 2? sin@ dédy (4.13) 


Using (4.9) in (4.13) 


aL, ce eee 
=> 20 (“sru?) | e 000?) ain Odo (4.14) 
0 
aL q 2 p2 2p2 
= 27 (“EoB0k?) ee x i ek Res? sin 6 dO (4.15) 
0 
Using the integral definition of modified bessel function (4.15) becomes 
a 2 
aL 2 p2 , AR 
=> In | —gBoR* | e* ™ x —— hh (KR)? 4.16 
( A gPo ) > ne (( ) ) ( ) 


Therefore we arrive at the expression for Total cross section 


((KR)) (4.17) 


Nile 


aL 2 , 
Oto = V2 (SE opar?) KRe&®) 7 


where Bo is the magnetic field strength of the solenoid, K is the incident photon 


momentum Ww is the frequency of the incident beam, L is the length of the solenoid and 
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R is the radius of the solenoid. 


The Axion transmittion probability Pirans can be found dividing o;,; by the lateral 


surface area of the solenoid of length L, and radius R 





(4.18) 


The QCD inspired axions with masses up to the ~ 1 eV range are our core motivation. 
We use the value of the coupling constant g from the CAST collaboration’s latest 
result to determine the magnitude of the entire cross-section. Following the coherent 
inverse Primakoff-effect, CAST is looking for axions formed in the sun and travelling 
to Earth by trying to detect photons from the conversion of axions inside a continuous 
magnetic field. CAST has set an upper bound on the magnitude of the axion-photon 
coupling constant of ggyy < 2.2 x 107"° GeV~! for an axion mass of mg S 0.4 eV, in 


~N 


collaboration with the Japanese axion helioscope Sumico . The (Axion Like Particles) 





10-6 
Exclusion domain for 
= pseudoscalar and scalar ALPs 
| 
> 
® 
O 
—™~ 10-7 
§ 
SH 
scalar 
10-8 
10-4 10-3 
Me / eV 


Figure 4.1. Exclusion limits for PS-ALP and S-ALP searches obtained in vacuum by 
the current OSQAR LSW experiment at 95% C.L. 


ALPS Collaboration’s exclusion limitations are also shown in Fig. 6, which summarises 
(Pseudo-Scalar-ALP) PS-ALPs and (Scalar-ALPs) S-ALPs searches. For the pseudo 
scalar and scalar searches for mg < 2 x 10~‘ eV, a limit for the diphoton couplings 
of Joyy < 3.5 x 107-§ GeV and gy, < 3.2x 107-8 GeV~'is found. These findings 
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impose the most restrictive limits on ALP searches in the almost massless limit that 
LSW experiments have yielded so far. Throughout this work, we will utilise the value 
P= G5, 0:88 107 GeV~'. Multiplying the 3D scattering cross section oj, with 
the incident photon flux F’ gives the number of events per second N;,,;. A few examples 


for the transmitted Axion conversion Probability, Cross section and Number of events 


are below in the Tables 4.1(a), (6), (c) 
















































































B [Tesla] | R[ cm] | oy: [barn] Nio:| sec” *] Ps 
10 1 2.699 x 107-7 | 9.906 x 107! | 4.296 x 1078° 
10 10 | 2.830 x 107? | 1.039 x 107!9 | 4.504 x 10782 
9 1.630 x 10~® | 5.982 x 107-14 | 1.297 x 10~-*4 
920 | 1.630 x 1072] 5.982 x 10719 | 1.297 x 1073! 
4.1(a) : For X-ray Photons (£ = 124 eV) 
B [Tesla] | R[ cm] | o¢¢ [barn] Nyot| sec] Piiss 
10 1 3.073 x 1074 | 1.128 x 107-2! | 4.891 x 107-42 
10 10 2.987 x 1078 | 1.096 x 10715 | 4.754 x 1073” 
2 7.075 x 107!8 | 2.596 x 10-2" | 5.630 x 10-*! 
20 6.324 x 1078 | 2.321 x 107!4 | 5.032 x 107% 
4.1(b) : For Visible light (E = 1.7 eV) 
B [Tesla] | R[ cm] | o¢¢ [barn] Neot| sec |] Pie 
10 1 8.702 x 10-27 | 3.194 x 107-84 | 1.385 x 1074 
10 10 |8.702 x 107?! | 3.194 x 10728 | 1.385 x 10-* 
2 2.005 x 10725 | 7.358 x 10788 | 1.596 x 107% 
20 | 2.005 x 107!9 | 7.358 x 10727 | 1.596 x 107% 





Table 4.1(a),(b).(c) : Total cross-section, number of events, and axion-photon con- 


version probability for various magnetic field strength (B), solenoid Radius (R) and 


4.1(c) : For Microwave Photon (E£ = 1.24 meV) 
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solenoid Length (LZ) = 10 cm values, as well as for g = 0.88 x 107!° GeV~t. The 
magnetic field units were converted from Tesla to eV using rationalised natural units, 
with the conversion being 1 T = 195 eV’. (please see appendix A in [63] for more 
details). 
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(a) B-R plot (b) K-B plot 





(c) K-R plot 
Figure 4.2 : Visualizing the behaviour of the total cross section o;.; with 
respect to arbitrary B, R, and K values 
(3D plot using matplotlib library in python - see Appendiz A for the detailed Graphs 
and Python codes ) 


We can see from the above 3D graphs that the arguments of the Modified Bessel func- 
tion Ja (( R)?) and the e~“*®)” are K and R. The behaviour of the total cross section 
is heavily influenced by these two arbitrary variables. The increase in solenoid radius 


R from 0 to 20 cm blows up the axion cross section by a factor of 7. Similarly, the 
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(a) K plot (b) R plot 


Figure 4.3 : Visualizing the behaviour of the total cross section o;,; with 
respect to arbitrary B, R, and K values 


similar pattern can be seen for the ranges of K and B. 


As a result, we’re going further into the equation (4.17), trying to figure out how 
the cross section and probability behave in relation to the variables K and R. In the 
absence of any change in R, we obtain the plot 4.3(a) by graphing the curve against the 
cross section. In the image, as one moves along the positive K axis, the slope steadily 
increases from 0 - 124 eV, indicating that with a constant solenoid radius, the rise in 
incident photon momentum accelerates the axion production rate represented by the 
increase in cross section. If we look at the slope of the R in figure 4.3 (b), we see that it 
is much larger than the slope of the K, which means that the increase in cross section 
with respect to an increase in solenoid radius is much greater than the increase in cross 


section with respect to incident photon momentum is relatively small. 


When you look at the equation 4.17, you will notice that the arguments of the expo- 
nential function and the modified bessel function are both (R* K)?, which results in a 
difference in the slope of the equation. As seen in figure 4.3(b), if we express frequency 
and momentum in energy units and simplify equation 4.17, we obtain the K? and R° 
terms, accordingly. The additional R? term contributes to the rise in the slope of the R 
VS Otor plot (figure 4.3(b)). If the solenoid’s radius is increased tenfold, it will result 
in an increase in the cross section oj on the order of eight times greater. When the 
incident photon moemntum is ten times greater, the cross section otot grows by an 
order of magnitude greater than two times. As a result, the slopes of the plots 4.3(a) 


and 4.3(b) deviate by such a wide margin from one another. 
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EXPERIMENTAL SEARCHES FOR 
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5.1 AXIONIC DARK MATTER 


Axions produced in the early universe with a mass ranging from eV to meV could 
account for the cold dark-matter component of the universe, as described in the review 
by P. Sikivie on cosmological axions in Sect. 1.6, axions produced in the early universe 
with a mass ranging from eV to meV could account for the cold dark-matter com- 
ponent of the universe axion detectors based on microwave resonant cavity studies are 
currently the most sensitive detectors for cold dark matter axions. The metallic con- 
ductor at the centre of such a device serves as a cavity for electromagnetic fields such 
as photons. It is possible to modify the resonant frequency of the magnetic-field per- 


meated cavity such that it coincides with the frequency of the axion field by adjusting 
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the size of the magnetic-field permeated cavity (which is related to the axion mass). 
Excess cavity power output could be detected with sensitive microwave receivers in 
this circumstance because the axions will resonantly change into quasi-monochromatic 
photons. Figure 5.3 depicts a schematic representation of the working principle of this 


type of detector. 


Axion distribution in the galactic halo would be represented by the width of the signal 
obtained, which would provide information not only about the axion mass but also 
about its distribution in the galactic halo. Furthermore, as illustrated in Fig. 5.3, the 
signal may possibly have finer substructure as a result of axions that have recently 
entered the galaxy and have not yet thermalized. The reader is directed to Sikivie’s 
review for a more extensive discussion of the axion distribution in the galactic halo, 


which is also included in this section. 


Consequently, because the dark-matter axion has an extremely small mass range, the 
expected signal is miniscule, and as a result the sensitivity of the experiment is critically 
dependent on the quality of the microwave receivers used in the experiment (see Chap. 
8 by Carosi and van Bibber). Small volume cavities (~ 1000 cm?) and high frequency 
transistor amplifiers were used in two pioneering experiments carried out in the 1980s at 
Rochester-Brookhaven-Fermilab (RBF) and the University of Florida (UF), but their 
sensitivities fell short by two or three orders of magnitude of the sensitivity required to 


probe the theoretically motivated region in the axion parametric spherical symmetry. 


In response to an increase in sensitivity, experiments such as ADMX were developed 
that are currently gathering data with a larger resonant cavity than the preceding 
studies, while still utilising the HFET amplifiers as previously described. The authors’ 
results are consistent with the absence of any axion signal, and they provide the best 
upper limits for the axion to photon coupling constant g,, in the lowest mass range 
(mg ¥ 10-*— 10~® eV). It is feasible to boost the sensitivity of ADMX even further 
by using a new amplifier technology based on Superconducting Quantum Interference 
Devices (SQUID), which is currently under development. This would allow researchers 


to probe the theoretically favored axion parameter range. 


The CARRACK experiment in Kyoto is another example of a second-generation ex- 
periment. Although it similarly makes use of a resonant cavity to generate the axion- 


converted photons, the signal detection is carried out using a Rydberg-atom single- 
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quantum detector, which has extremely low noise photon counting capabilities due to 


its low noise. This experiment is still in the early stages of development. 


A suitably light axion, whose density relative to the Universe’s critical density is given 


by, 
6 
Qa x (=) (5.1) 


is a promising dark matter candidate. The whole dark matter density of the Universe, 





Qn * 0.27, would be accounted for by an axion of mg * 20 eV (within a factor 
of nearly two). Much lighter axions would overclose the Universe if the initial mis- 
alignment angle was not tuned; hence, mg ~ 1 eV may be regarded a strong lower 
limit on the axion mass mg. The relative contribution of axion creation via the vac- 
uum realignment mechanism and radiation from topological flaws has been a source 
of debate (axion strings, domain walls). Although a definite answer is doubtful in the 
near future, the prevailing agreement is that the sum of all contributions to Q, might 
raise the axion mass corresponding to Qpyy by two orders of magnitude. Because we’re 
mainly concerned in creating conservative ranges for experiments right now, searches 
should aim for the eV scale. The formation of thermal axions in cosmology provides 
an upper bound on axion mass, the hot dark matter limit of around 1 eV, however we 


are not interested in this topic. 


Axions are likewise subject to stringent restrictions as a result of stellar evolution. Ax- 
ions with masses greater than ~ 16 meV would have quenched the SN1987a neutrino 


pulse, therefore constraining the axion mass from above. 


Galactic globular clusters provide the best indirect stellar bound on the axion-photon 
coupling, g¢yy < 0.66 x 10-'°GeV~!. This bound just eclipsed the greatest direct 
stellar bound, gg, * 0.88 x 10-!°GeV~ , obtained by the CERN Axion Solar Telescope 
(CAST) hunt for solar axions. Such couplings for PQ axions correspond to masses well 
outside the open mass region (i.e., 107° eV < mg < 10-7 eV). A brief literature 
has recently indicated that the brightness function of white dwarfs (degenerate stars 
undergoing gravothermal cooling) may require an extra cooling process, which might be 
explained by axions in the ~ 10 meV range. The evidence is inconclusive; additionally, 


such masses are beginning to encroach on the SN1987a favored zone. 


The open mass window for axions has long been thought to be in the 10~® — 10? eV 
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Figure 5.1. Ranges of exclusion as described in the text. The dark gaps represent the 
approximate CAST and ADMxX search ranges, while the green sections represent the 
future upgrade’s intended reach. Using z = 0.56 and the KSVZ values for the 
coupling strengths, coupling strength limitations are translated into limits on mg and 
fe . For regular or modified axions, the ” Beam Dump” bar is an approximate 
approximation of the exclusion range. The DFSZ model with an axion-electron 
coupling of cos? 8 = 1/2. is used in the ” Globular Clusters” and ” White Dwarfs” 
ranges. The exclusion range for Cold Dark Matter is particularly ambiguous; ranges 
for pre-inflation and post-inflation Peccei-Quinn transitions are presented. 


range. String theories, on the other hand, are abundant in axions or ALPs greater than 
100 in any given realisation yet such theories naturally favour fy ~ 10'°— 10'® GeV, 
corresponding to neV scale masses. It’s impossible to identify which of them, if any, 


solves the strong-CP problem and which is cosmologically relevant. 


5.2 AXION HELIOSCOPES 


Axion helioscopes, or magnetic solar telescopes, are currently the most sensitive axion 
experiments in the mass range of 10-°eV < mg < leV. The underlying physical basis 
of such a telescope is represented in Fig. 5.2 and is based on an idea proposed by 
P. Sikivie in 1983. If axions are created in the Sun, they would arrive on Earth as a 
virtually parallel axion beam after about 500s (the apparent diameter of the axion 
emission zone on the solar disc is ~0.1°). The time reversed Primakoff effect, i-e., 


a+ Wirtual — Y, is used to detect the axion on Earth, where the axion interacts with 
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Figure 5.2: An axion helioscope’s operation is based on the following concept. In a 
high magnetic field, axions generated in the Sun’s core by the Primakoff effect would 
be transformed back into photons, completing the cycle. These photons, which have 
an energy spectrum that matches that of the arriving axions, may be captured by an 

X-ray detector located at the end of the magnet field area in due course. 


a virtual photon provided by a transversal magnetic field and reconverts into a real 
photon. Because these ” reconverted” photons have the same energy and momentum as 
the axion, their energy distribution is identical to the solar axion’s energy spectrum. An 
X-ray detector on the magnetic field’s end would gather photons and look for an axion 
signal above the detector’s background. For the given differential solar axion spectrum, 
the estimated number of reconverted photons reaching the detector is For the given 
differential solar axion spectrum, the estimated number of reconverted photons reachng 
the detector is. 
d®, 
NN = / ante Stdk (5.2) 


The differential axion flux at the Earth is d®,/d£, the effective axion-sensitive mag- 
netic aperture is S, the measurement period is t, and the conversion probability of an 
axion into a photon is P,_,y, the magnet is positioned in the same way as traditional 


telescopes, allowing it to monitor the Sun’s path across the sky. 


5.2.1 Pioneering Axion Helioscope 


A 1.8 m long magnet driven by a 2.2 T transverse magnetic field was used to build the 
first axion helioscope in the early 1990s by the RBF collaboration . The photons from 


the axion-to-photon conversion were collected using a proportional counter running 
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with a gas mixture of argon and methane (90% to 10%) as a detector. For a period of 
~ 15 min day ~', the helioscope may be directed to the Sun. This configuration was 
used to investigate two distinct areas of the axion mass range (with and without a gas 


in the conversion volume), placing the following constraints on gy, and mg 


Géyy < 3.6 X 10°? GeV" for mg < 0.03 eV 
Gbyy < 7.7 X 107? GeV" for 0.03 eV < mg < 0.11 eV 


5.2.2 The Tokyo Axion Helioscope (Sumico) 


In the late 1990s, Tokyo constructed a new helioscope with enhanced sensitivity (2.3 m). 
For this experiment, sixteen PIN photodiodes were employed as X-ray detectors, and 
the superconducting magnet was mounted on a polar mount, allowing the system to 
track the Sun in a declination range of —28° to +28° and a right ascension range of 360°. 
The Tokyo helioscope was able to watch the Sun for 24 h a day thanks to a superior 
tracking mechanism, and the Sun observation period was substantially enhanced. The 
Tokyo cooperation was able to employ roughly 50% of their instrument’s duty cycle to 
follow the Sun, with the remainder of the time being used for background observations. 
The experiment was separated into two parts, similar to the original method of the 
first phase was completed in 1997 with an evacuated conversion zone, and the second 
phase was completed in 2000 with the magnet pipe filled with helium gas to increase 
the sensitive axion mass range. The top limits on g¢,y derived from both data taking 
periods are 


Jovy < 6X 107'’GeV-! for mg < 0.03eV 
Jovy < 6.8 — 10.9 x 10°"°GeV~" for mg < 0.3eV 


5.2.3. CERN Axion Solar Telescope (CAST) 


The CERN Axion Solar Telescope, the world’s most sensitive axion helioscope, is now in 
operation at CERN (CAST). It was constructed in 2002 by refurbishing a 9.26 m long 
prototype of a LHC dipole magnet with dual apertures. A transverse dipole magnetic 
field is homogeneously interleaved between the two parallel pipes of the magnet, each 


with a bore area of A = 14.5 cm?. The highest magnetic field that this prototype 
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magnet can produce is ¥ 9 T. During dawn and sunset, the magnet is positioned on 
an azimuthal movable platform that allows it to track the Sun for 1.5 h(+8° in height 
and = 100°in azimuth). Three detectors are attached to the ends of the magnet pipes 
to detect axions: a gaseous detector with a novel MICROMEGAS readout , an X-ray 
telescope with a pn-CCD chip as focal plane detector to detect axions during sunrise, 
and a Time Projection Chamber (TPC) that covers the opposite end of the magnet 
pipes, ready to collect axions during sunset. During the years 2013 to 2015, CAST 
collected data with the conversion region removed. An upper limit on the axion to 


photon coupling constant has been determined in the absence of any axion signal. 
Joyy < 0.66 x 107'° GeV" for mg < 0.02 eV 


at the 95% confidence level. This is the best laboratory limit a helioscope has ever 
attained. It encompasses a wide range of masses and outperforms the best astrophys- 
ical limit, which was determined from the development of horizontal-branch stars in 
globular clusters . The CAST experiment’s design was updated in 2005, allowing it to 
proceed to the second phase of the experiment, when the magnet pipes are filled with 
“ He or He gas. CAST will be able to increase its sensitivity to a maximum axion 


mass of © 1.1 eV in this setup. 


5.2.4 Orbiting X-ray Telescopes 


As with the axion helioscopes discussed above, the same operating concept applies to 
orbiting detectors sensitive to hard X-rays (for example, for the RHESSI solar X-ray 
observatory, which has a threshold greater than ~ 3 keV). Axion-to-photon conversion 
can occur either in the terrestrial magnetic field or in the magnetic field near the solar 
atmosphere, depending on the circumstances. Despite the significant disparities be- 
tween the two methods, their sensitivity can compete with the finest earth-bound helio- 
scope (see the next section), but only for axion rest masses in the range of 10~4eV /c? or 
less. According to the findings of this article, solar X-ray telescopes in space may be 
used as the most sensitive solar axion antennas for large axions (mg > 10~* eV/c?) by 
employing the solar magnetic fields near the photosphere in addition to the solar X-ray 
telescopes. One may be tempted to think of such a system as an indirect approach in 


axion helioscopy, which would be incorrect. When compared to direct axion-detection 
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approaches utilising axion helioscopes, the findings acquired in this manner can be just 


as significant, if not more so, and eventually have a higher built-in sensitivity. 


Table 1.Axion helioscopes are being compared in terms of their figures of merit . When 
calculating the magnetic field strength, B equals the strength of the magnetic field, L 
equals the length, S equals the effective axion-sensitive magnetic aperture, and ¢ is the 


tracking time per day (for the orbiting telescopes). 








Helioscope (BL)? (BL)*S (BL)*St 

Tt? T? m+ T? m‘ hours 
RBF 16 ~3x 107? ~1x 107? 
Sumico 85 —10x10-2 ~120x 10? 
CAST 6946 2000 x 10-2 6000 x 107? 
In orbit 324 20000 x 10-2 — 





5.3 AXION MICROCAVITY EXPERIMENT 


In the solar vicinity, the black halo of our Milky Way galaxy has a density of order 
10~*4 er/cm?. The halo particles have v velocities in the 10~c range. We are sur- 
rounded by a pseudo-scalar field pulsating with angular frequency if the dark matter 
is axions: 

gS Bg Wigs sng” = mg (1+ 0 (10°*)) (5.3) 


The axion electromagnetic interaction (4) becomes 


La = -00 ak - Bo (5.4) 
in the presence of an externally supplied magnetic field Bo. As mentioned in the 
preceding section, it enables for the conversion of axions to photons and vice versa. It 
is advantageous to have the conversion process occur inside an electromagnetic cavity 
in the case of dark matter axions, assuming their mass is in the 10~° to 10~4 eV range. 
When one of the cavity modes equals the axion signal’s angular frequency, the cavity 


collects the photons produced and boosts the conversion process through resonance. 
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The conversion power from an optimised experiment is given by 
Pi VB i Oye a 
a ~Ioyy 00a I nun (On, Qa) 
V Be gy \2 
=0.5 x 10° w(=— | (=) c(42) 
. (sr) (= | 0.36 
S Qa 
(a x 10-24 ¢ =) 


) min (Qs. Qa) 


(5.5) 





x ( ms 

27 (GHz) 
V denotes the cavity volume, By denotes the magnetic field, Q; denotes the cavity’s 
loaded quality factor (defined as centre frequency over frequency bandwidth), Qa = 
10° denotes the quality factor of the axion signal (axion energy over spread in energy 
or 1/67), 0q denotes the axion mass density at the detection point (earth), where 
Cimn is the form factor for one of the transverse magnetic cavity modes (TMimn). 
Essentially, this form factor is the normalised overlap integral of the external static 
magnetic field, Bo(x), and the oscillating electric field, E,,(x)e“", of that particular 


cavity mode. It can be determined with the use of 


L. | Sy GrE,, - Bol 
«BBV f, d®ze|E.|? 





(5.6) 


in where € denotes the dielectric constant of the cavity. The TMoio mode has the 
highest form factor (Co19 = 0.69[19]) for a cylindrical cavity with a homogenous lon- 
gitudinal magnetic field. Though Model-dependent, equation (5.5) can offer an idea of 
the exceedingly small signal expected from axion-photon conversions in a resonant cav- 
ity, measured in yoctowatts (10774 W). This is far less than the 2.5 x 107?! W of power 
received from the Pioneer 10 spacecraft’s 7.5 W transmitter’s last transmission in 2002, 
when it was 12.1 billion kilometres from Earth. Furthermore, because the experiment 
would require tuning orders of magnitude of frequency in short frequency steps, the 
Dicke radiometer equation governs how long one can integrate at each frequency to 


optimise signal-to-noise: 
P, Ps Avt 
22a (5.7) 
n Py kT, 


The signal-to-noise ratio is s/n, the signal’s bandwidth is Av, the integration time is 





t, and the signal and noise powers are P;and P,, respectively. Ty = Tpnys + Tetec - 
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Figure 5.3. The microwave cavity search for dark-matter axions is depicted in this 
diagram. In a strong magnetic field, axions resonantly convert to a 
quasi-monochromatic microwave signal; the signal is collected from the cavity by an 
antenna, amplified, down converted, and the power spectrum estimated by an FFT. 
Any fine structure in the signal would provide crucial details about our galaxy’s 
formation. 


ADMxX recently completed a transition from conventional heterojunction field-effect 
transistors (HFETs or HEMTs) with a noise equivalent temperature of Tyee ~ 2 K, to 
superconducting quantum interference device (SQUID) amplifiers, whose noise equiva- 
lent noise temperatures can reach the quantum limit, Tiree ~ 50mK at 750 MHz when 
cooled to comparable physical temperature limit. This method will allow them to reach 
the DFSZ model axions and cross the open mass range considerably more quickly in 


the long run. Carosi describes his microwave cavity studies in this issue of Focus. 


5.3.1 Early experiments 


Soon after the experimental concept was published, two pilot initiatives in the (1 — 
4 GHz) range were launched at Brookhaven National Laboratory (BNL) and the Uni- 
versity of Florida (UF). These used a few litre cavities and the greatest conventional 
amplifiers available at the time, such as HFET amplifiers. These amplifiers were not 
sensitive enough to reach PQ axions with noise temperatures in the ~ 3 — 20K range. 
Nonetheless, these two experiments laid the groundwork for much of the design philos- 


ophy and knowledge regarding microwave cavities that is still used in today’s research. 
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In the history of microwave cavity experiments, the CARRACK experiment in Kyoto, 
Japan, was a watershed moment. The goal of CARRACK was to significantly lower 
the system noise temperature by lowering the physical temperature of the experiment 
to ~ 15mK with a *He — “He dilution refrigerator and eliminating the amplifier noise 
contribution by using a Rydberg-atom single-quantum detector instead of a standard 
linear amplifier. The standard quantum limit (SQL) is an irreducible noise contribution 
that occurs in linear amplifiers: kTgq, = hy. The Rydberg-atom single-quantum 
detector can be thought of as a tunable radio-frequency (RF) photomultiplier tube, 
with the photon interacting as a particle rather than a wave, avoiding the SQL. Tada 
et al. observed the cavity’s blackbody photon spectrum as a function of temperature 
from 2.527 GHz to T = 67mK, about a factor of two lower than the SQL of ~ 120mK. 
CARRACK was a success from a technological standpoint, but it was eventually too 


complicated to be used as a production experiment. 


5.3.2 ADMX- Axion Dark Matter eXperiment 


The Axion Dark Matter eXperiment (ADMX) at Lawrence Livermore National Labo- 
ratory is seen in Fig. 8.1 as a schematic representation. It comprises of a cylindrical 
copper-plated steel cavity with two axial tuning rods that is housed within the cavity 
of the experiment. As a result, the resonant frequency may be perturbed by moving 
them transversely from the border of the hollow wall to its centre. The cavity itself is 
housed within the bore of a superconducting solenoid, which generates a powerful and 
continuous axial magnetic field around the cavity itself. In this case, a tiny adjustable 
antenna is used to link the cavity’s electromagnetic field to the low-noise reception 
electronics. First, two cryogenic amplifiers with extremely low noise are connected in 
series to amplify the signal in this circuit. Using a room-temperature post-amplifier, 
the signal is then amplified once more before being sent to a double-heterodyne re- 
ceiver. It comprises of an image reject mixer that reduces the signal frequency from 
the cavity resonance (hundreds of MHz — GHz) to an intermediate frequency (IF) 
of 10.7 MHz, which is then amplified. In order to reject noise power outside of a 
35 kHz window centred on the IF, a crystal bandpass filter is utilised. Once this 
has been accomplished, the signal is blended down to practically audible frequencies 
(35 kHz) and examined by fast-Fourier transform (FFT) electrical circuitry. The FFT 
circuitry computes a 50 kHz bandwidth that is centred on 35 kHz. By adjusting the 
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(variable coupling) 


Microwave cavity 


Figure 5.4.The ADMX experiment is depicted in a schematic design that includes 
both the resonant cavity (which is housed within the bore of a superconducting 
solenoid) and the receiver electronics chain. 


tuning rods to produce a new resonant TM 919 mode, data is collected about every 
1 kHz. Some of these components will be discussed in further detail in the following 


sections. 


The primary magnet for ADMX was created in such a way that the B3V contribution 
to the signal power was maximised (Figure 5.5). In the end, it was concluded that a 
superconducting solenoid would be the most cost-effective option, and that its extraor- 
dinarily large inductance (535 Henry) would have the additional benefit of maintaining 
a very steady magnetic field. An open magnet bore in the cryostat (see Fig. 5.5) allows 
the experimental insert with the cavity and liquid helium (LHe) reservoir to be lowered 
into the 6t magnet coil. The magnet coil is housed in a 3.6 m tall cryostat (see Fig. 
5.5) with an open magnet bore. During operations, the magnet itself is immersed in 
a 4.2 K LHe bath in order to maintain the superconductivity of the niobium-titanium 
windings on the magnet. For the most part, the magnet was held at a field strength 
of7.6 T in the solenoid centre (with a drop to roughly 70% at the ends), although the 
experiment has recently been conducted at field strengths of up to 8.2 T Because of the 
utilisation of cylindrical cavities in the ADMX experiment, it is possible to optimise 


the axion conversion volume in the solenoid bore. A copper-plated steel cylinder with 
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Figure 5.5: (a) ADMX hardware overview, including the superconducting magnet 
and cavity insert (b) Superconducting solenoid “bucking” magnet. 





sealed ends serves as the basis for these cylinders. The Helmholtz equation, 

V7O+ hb =0 (5.8) 
where the wavenumber k is given by 

k? = pew — 6? (5.9) 


where § is the eigenvalue for the transverse (x,y) component, can be used to deter- 
mine the electromagnetic field structure inside a cavity. The cavity modes are the 
standing wave solutions to (5.8). The boundary conditions of an empty cavity allow 
only transverse magnetic (TM) modes (B, = 0) and transverse electric (TE) modes 
(E,=0) Three-dimensional standing waves are represented by the TMjmnn modes, 
where | = 0,1,2,... is the number of azimuthal nodes, m = 1,2,3,... is the number 
of radial nodes, and n = 0,1,2,... is the number of axial nodes. The axions have the 


strongest coupling to the TMo19 mode of the simplest structure. 


By introducing metallic or dielectric tuning rods placed axially into the cavity, the 
resonance frequency of the T’Moig mode may be changed. As you move closer to the 


centre, metallic rods boost the cavity resonant frequency, whereas dielectric rods reduce 
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it. These rods are connected to the ends of alumina arms of ADMX, which pivot around 
axles located in the upper and lower end plates (Fig. 5.5). The tuning rods are swung 
in a circular arc from the cavity border to the centre by stepper motors installed at 
the top of the equipment (see Fig. 5.4). The stepper motors are connected to a gear 
reduction that converts one step into a 0.15 arcsecond rotation, which corresponds to 
a shift of~ 1kHz at 800MHz resonant frequency. 


800 





f (MHz) 
Frequency 





Tuning rod steps 


3 3.5 4 45 5 5.5 6 
Tuning rod steps (x 10°) 


Figure 5.6 .A cavity with two copper tuning rods has a mode structure. Left : 
Resonant mode frequencies as measured with a swept RF signal when one tuning rod 
is held at the cavity edge and the other is moved toward the centre. Right : sketch 
of a mode crossing. 


TEM modes (B, = FE, = 0) can be supported in the cavity by the inclusion of metallic 
tuning rods. They do not couple to the axions, but they can weakly couple to the 
vertically mounted receiver antenna, just as the TE modes (due to imperfections in 
geometry, etc). The numerous resonant modes alter as a copper tuning rod is moved 
from near the cavity wall to the centre, as shown in Figure 5.6. The tuning-rod location 
has no effect on the TEM and TE modes, while the TM modes increase in frequency 
as one of the copper rods approaches the cavity centre. As a result, areas where a 
TM mode crosses a TE or TEM mode appear (referred to as mode mixing). These 
mode mixings (seen in the right section of Fig. 5.6) bring in new information Gaps in 
frequency that can’t be scanned. To solve this, the cavity was eventually filled with 
LHe, which lowered the microwave index of refraction to 1.027, decreasing the mode 
crossings by 2.7% and permitting scanning of previously unaccessible frequencies. The 
ADMxX experiment will be upgraded with a dry system in which both tuning rods are 


adjusted to ”dance” around the mode crossover without the need of LHe. 
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The quality factor Q, which is a measure of the sharpness of the cavity response to 
external excitations, is a crucial characteristic of the resonant microwave cavity. It is 
a dimensionless variable that may be expressed in a variety of ways, such as the ratio 
of stored energy U to power loss P;, each cycle: Q = woU/P,. Sweeping a RF signal 
through the weakly-coupled antenna in the cavity top plate determines the quality 
factor Q of the TMo1o mode (see Fig. 5.4). In most cases, the cavity’s unloaded Q 
is ~ 2 x 10°, which is extremely close to the theoretical maximum for oxygen-free 
annealed copper at cryogenic temperatures. The insertion depth of the main antenna 
is modified during data collection to ensure that it fits the cavity’s 50 Q impedance 
(called critically coupling). Half of the microwave power in the cavity reaches the 
electronics via the antenna when the antenna is critically linked, while the other half 
is dissipated in the cavity walls. Overcoupling the cavity lowers Q and therefore limits 
signal amplification, whereas undercoupling the cavity limits microwave power entering 


the electronics. 


5.3.2.1 Gen 2/Phase Ila upgrades 


The Gen 2 ADMX experiment will look for axion dark matter down to DFSZ coupling 
in the mass range of 2 — 40 preV (0.5 — 10 GHz). Several multiyear RD projects have 
been formed to achieve this goal. The Cavity Working Group is working on designing 
and fabricating next-generation microwave cavities for higher-mass searches over the 
long term. To lower heat loads, piezoelectric motor drives are being tried, and tunable 
high-frequency SQUID amplifiers and Josephson parametric amplifiers (JPA) are being 
created. High-field magnets are also being studied and costed as part of the effort. 


A tiered frequency mechanism will be used to reach search capabilities of more than 
one order of magnitude. The frequency ranges include 0.5 — 1 GHz, 1 — 2 GHz, 2 — 
4 GHz,4—6 GHz,6—8 GHz, and 8—10 GHz. The present ADMX cavity will cover the 
first frequency range. The Pound stabiliser system will be used to frequency lock four 
cavities in the second frequency tier. The frequency locked cavities used in the 2-4 GHz 
and 4-6 GHz tiers will be ~ 16 and ~ 32, respectively. A photonic band-gap cavity, a 
complicated system of many tuning posts that move in unison to maintain symmetry 
to enhance coupling, C’, will be required for the 6 — 8 GHz range.The research on the 
ultimate tier configuration is still in its early phases. To achieve DFSZ sensitivity at 


greater masses, complex manufacturing techniques are used to eliminate mode mix- 
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Figure 5.7 .a) A graphic representation of the cavity design for 1-2 GHz. The cavities 
will be frequency locked using the Pound stabiliser. b) At the top, a shot of a 
photonic band-gap prototype being built. Bottom: the electric field strength as a 
result of a simulation. The electric field will be significantly localised as symmetry is 
broken, resulting in a drop in C. c) Prototype of a piezoelectric motor drive. 


ing and symmetry breaking, two abnormalities that drastically diminish C’ within the 
search range, while also increasing the quality factor of the cavities. Hybrid supercon- 
ducting cavities are also being researched. Hybrid cavities would use superconducting 
thin-film coverings on the vertical surfaces to boost the detector’s Q by a factor of ten. 
In addition, to improve scan reliability, in situ mode identification techniques and field 
characterization approaches are being tried. The 4-cavity concept is depicted in Figure 
5.7(a — b), which includes a graphic image of a photonic band-gap cavity as well as a 
photograph taken inside a prototype. To address the experiment’s fine motion control 
requirements, ADMX currently employs worm gears. The gears generate a lot of heat, 
which will be too much for a large number of cavities. In order to prepare for future 
searches, piezoelectric motor drives are being tested on the sidecar. A snapshot of the 


motor drive prototype is shown in Figure 5.7(c). 


5.3.3 Future of Microwave Cavity Experiments 


Several new search concepts for detecting dark matter axions have emerged in recent 
years. The NMR-based CASPEr was created to investigate the mass range 10~° — 


10-°eV. Other potential ideas, on the other hand, deserve a quick notice. 


Rybka et al, have reported an experiment based on an open microwave resonator 
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structure suitable for axions in the 100 — 1,000 peV decade, with a microwave pho- 
ton wavelength ranging from 1 to 10 mm after conversion. Because the conversion 
probability is much increased when the applied magnetic field B(r) follows the photon 
mode’s electric field E'(r), a wiggler-like magnetic field with a continuously controllable 
periodicity is required to sweep out the mass range. Rybka et al. achieved this by em- 
ploying a sequence of wire planes in an open resonator of alternating sign in current; 
these planes will eventually be made up of superconducting wires or stripes imprinted 
on a thin substrate. With couplings > 4 x 10-’ GeV~!, the early prototype was able 
to rule out axions as dark matter between 68.2 and 76.5 peV. Finally, below DFSZ 


axions for that mass range, sensitivity to couplings of 10~!° GeV~! should be possible. 


Sikivie et al, have focused on the opposite issue of the microwave cavity experiment, 
namely extending the search lower in mass, where the cavity and the magnet that 
encloses it become unfeasibly huge below 107° eV(~ 250 MHz). The authors rec- 
ommended replacing the cavity with a lumped-parameter LC circuit external to the 
magnetic field, stimulated by a pickup loop inside the magnet threaded by the mode’s 
transverse magnetic field, to get around this constraint. According to the authors, 
this technique would be most sensitive around 10~" eV, probing the band of PQ axion 
models, with a magnet similar to that of the existing ADMX. Such a lumped parameter 
LC circuit search for dark matter axions, if fully conceived and executed, might bridge 
a major gap where the cavity-based experiment and the NMR-based experiment may 


have difficulties overlapping. 


The proposal for the Cosmic Axion Spin Precession Experiment (CASPEr) intends to 
pave the way for a new route in the experimental hunt for axion dark matter. CASPEr 
will employ nuclear magnetic resonance (NMR) methods to identify the spin precession 
generated by axion dark matter. This unique technique adds to prior efforts: ADMX is 
sensitive to higher axion frequencies, CASPEr, on the other hand, will cover the lower 
frequencies where the axion emerges from energies f, ~ 10'° GeV — 10!° GeV This 
range is extremely difficult to approach with any other technique, while certain astro- 
physics techniques may be able to investigate it . A detection in such an experiment 
would not only mark the discovery of dark matter, but would also reveal insights into 
the high-energy scales from which such an axion might emerge, close to fundamental 


scales like the grand unification, Planck, or string scales. 
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CHAPTER 6 


SUMMARY & CONCLUSION 


This article conducts an in-depth examination of the realm of axionic dark matter. 
After addressing a brief history of axions, their cosmological origins, existing models, 
and the rationale for considering axions as cold dark matter in the first two chapters of 
this article, we were able to focus on studying the electromagnetic coupling of axions via 
the primakoff effect and calculating the three-dimensional cross-section for the scattered 
axion transmitted out of a Gaussian distributed magnetic field centred around the axis 
of a solenoid. These findings are expanded upon and used to create a comparative 
study of how fine tuning different parameters such as the magnetic field (B) and the 
radius of the solenoid (R) affects the cross section, number of events, and transmission 


probability of scattered axions. The results are presented in a tabular format in Table 
4.1 (a), (b) ,(c). 


These findings indicate that raising the energy of the incident photon from microwave 
to x-ray,by 5 orders of magnitude, also increases the transmission probability by ap- 
proximately 19 to 17 orders of magnitude, depending on the strength of the magnetic 
field (B). As the energy level increases, the correlation between Raidus and the ax- 
ion transmission probability Pian; weakens. i,e A tenfold increase in the radius R of 
the solenoids raises the Ptrans by five orders of magnitude for visible and microwave 


photons, but only by three orders of magnitude for X rays. 


By examining the behaviour of equation 4.17, which is the transmitted axion cross- 
section oj, using Python data analysis and graphing libraries such as matplotlib, 
numpy, and scipy, we conclude that a tenfold increase in the K,R, and B values results 
in a sevenfold increase in the cross-section. This is due primarily to the modified bessel 
function and the exponential term with K and R arguments. Further decomposing 
equation 4.17 by holding all other variables constant and adjusting only the K and 
R variables one at a time and plotting it against the cross section oj; reveals that 
the slope difference between both of these plots presented in figure 4.3 allows us to 


fine-tune the experiment for the desired axion conversion probability. 
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Chapter 5 presents a full review of the axion studies based on the inverse primakoff 
principle. We address the quest to discover the solar axion through a variety of studies 
and collaborations. Using axion telescopes such as CAST and SUMICO, as well as X- 
ray telescopes, which investigates the limitations of the axion photon coupling constant 
and mass. CAST recently released a work in which they established an upper limit of 
Jovy < 0.66 x 1071? GeV_, for mg smaller than 0.02 eV with a 95 % confidence level. 


This is the highest level of performance yet achieved by a laboratory limit helioscope. 


Additionally, this thesis discussed cutting-edge research being conducted on Axion 
Microcavity Experiments such as BNL, CARRACK, and ADMX, as well as their up- 
coming enhancements. Microcavity experiments may now explore a range of 10° to 10° 
eV using current technology. ADMX and CASPEr will work in tandem to cover higher 
and lower accident frequencies that are difficult to approach with other techniques. 
This will increase the likelihood of Axion discovery in the near future and may also aid 


in probing more fundamental energy scales. 


It’s important to remember that finding the axion would not only answer the Strong- 
CP problem and shed light on the nature of dark matter, but it might also open up new 
vistas in astrophysics, cosmology, and quantum physics. There may be fresh insights 
into the creation of the Milky Way if fine structure is discovered in the axion spectrum. 
In fact, because of the axion’s large de Broglie wavelength (A, ~ 10-100 m), fascinating 
quantum experiments might be carried out at macroscopic scales . The axion is a very 
intriguing dark matter candidate because of all of these intriguing possibilities that are 


within reach of present and near-future technology. 
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APPENDIX A 


Python Plots and Codes 


A.1 3D plots 


A.1.1 oj vs B vs R plot 


PYTHON CODE : 








#importing the essential libraries 
from mpl_toolkits import mplot3d 
import scipy.special as spl 

import numpy as np 

import matplotlib.pyplot as plt 


import matplotlib.ticker as mtick 


#assigning the constant values for the cross section plot; 
w=124 

L=10 

g=8.8E-20 

y=124 


#defining our cross section function as z(x,y); 
def z_function(x, B): 
return (np.sqrt(2*np.pi)*(Cy*L*g*B* (pow(x,2)))**2 * (1/16))*y*pow() 
x, 5)*np.exp(-(y**2) *(x**2))* 
spl.iv(0.5, Cy ¥**2)*(x**2))) 


#defining our variables- R=x and B=y; 
x = np.linspace(0O, 0.2, 100) 
B = np.linspace(0, 1950, 100) 


#mapping the function z(x,y) to meshgrid Z(X,Y) 


X, Y = np.meshgrid(x, B) 
Z = z_function(X, Y) 
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A.1. 3D plots 


#plotting wireframe plot Z(X,Y) 

fig = plt.figure(figsize=(20,20)) 

ax = plt.axes(projection="3d") 

ax.plot_wireframe(X, Y, Z, cmap='autumn') 

ax.view_init (32.5,-75) 

#labeling the axis for wireframe 

ax.set_xlabel('solenoid Radius (R)'+ r'$\rightarrow$') 
ax.set_ylabel('Magnetic field (B)'+ r'$\rightarrow$') 
ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#fixing the figure size and zaxis format 
fig = plt.figure(figsize=(20,20)) 
ax = plt.axes(projection='3d') 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#plotting Z(X,Y) with color bar 

surf=ax.plot_surface(X, Y, Z, rstride=1, cstride=1, 
cmap='rainbow', edgecolor='none') 

fig.colorbar (surf , shrink=0.5,aspect=7) 

ax.view_init (32.5,-75) 

#labeling the axis for surface plot 

ax.set_xlabel('solenoid Radius (R)'+ r'$\rightarrow$') 

ax.set_ylabel('Magnetic field (B)'+ r'$\rightarrow$') 

ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


#fixing the figure size and zaxis format 
fig = plt.figure(figsize=(20,20)) 
ax = plt.axes(projection="3d") 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#plotting Z(X,Y) with color bar 

qcs=plt.contourf(X, Y, Z,20, cmap='rainbow') 
plt.colorbar(qcs,shrink=0.5,aspect=7) 

ax.view_init (32.5,-75) 

#labeling the axis for countour plot 

ax.set_xlabel('solenoid Radius (R)'+ r'$\rightarrow$') 
ax.set_ylabel('Magnetic field (B)'+ r'$\rightarrow$') 
ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 
#display our 3d plots 

plt.show() 
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A.1. 3D plots 


OUTPUT : 





Figure 7.1(a) : 
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Figure 7.1(b) 


58 


0200 


wireframe. 





0175 
0.200 


surface. 


Fo) 
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Figure 7.1(c) 
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A.1. 3D plots 


A.1.2 oj vs K vs R plot 


PYTHON CODE : 








#importing the essential libraries 
from mpl_toolkits import mplot3d 
import scipy.special as spl 

import numpy as np 

import matplotlib.pyplot as plt 


import matplotlib.ticker as mtick 


#assigning the constant values for the cross section plot; 
B = 1950 

w=124 

L=10 

g=8.8E-20 


#defining our cross section function as z(x,y); 
def z_function(x, jy): 
return (np.sqrt(2*np. pi) *(Cy*L*g*Bx* (pow (x,2))) **2*(1/16)) *y*pow(x, 
5) *np.exp(-(Cy**2) *(x**2))*spl. 
iv(0.5, (Cy **2) *(x**2))) 


#defining our variables- R=x and K=y; 
x = np.linspace(0O, 0.2, 100) 
y = np.linspace(0, 125, 100) 


#mapping the function z(x,y) to meshgrid Z(X,Y) 
X, Y = np.meshgrid(x, y) 
Z = z_function(X, Y) 


#plotting wireframe plot Z(X,Y) 

fig = plt.figure(figsize=(20,20)) 

ax = plt.axes(projection="3d") 

ax.plot_wireframe(X, Y, Z, cmap='autumn') 

#labeling the axis for wireframe 

ax.set_xlabel('solenoid Radius (R)'+ r'$\rightarrow$') 
ax.set_ylabel('Photon momentum (K)'+ r'$\rightarrow$') 
ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#fixing the figure size and zaxis format 


60 








A.1. 3D plots 


fig = plt.figure(figsize=(20,20)) 

ax = plt.axes(projection='3d') 

plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 

ax.view_init (32.5,-75) 

#plotting Z(X,Y) with color bar 

surf=ax.plot_surface(X, Y, Z, rstride=1, cstride=1, 
cmap='rainbow', edgecolor='none') 

fig.colorbar (surf ,shrink=0.5,aspect=7) 

ax.view_init (32.5,-75) 

#labeling the axis for surface plot 

ax.set_xlabel('solenoid Radius (R)'+ r'$\rightarrow$') 

ax.set_ylabel('Photon momentum (K)'+ r'$\rightarrow$') 

ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


#fixing the figure size and zaxis format 
fig = plt.figure(figsize=(20,20)) 
ax = plt.axes(projection="3d") 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#plotting Z(X,Y) with color bar 

qcs=plt.contourf(X, Y, Z,20, cmap='rainbow') 
plt.colorbar(qcs,shrink=0.5,aspect=7) 

ax.view_init (32.5,-75) 

#labeling the axis for countour plot 

ax.set_xlabel('solenoid Radius (R)'+ r'$\rightarrow$') 
ax.set_ylabel('Photon momentum (K)'+ r'$\rightarrow$') 
ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


#display the 3d plots 
plt.show() 
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A.1. 3D plots 


OUTPUT : 
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A.1. 3D plots 


A.1.3 oi vs K vs B plot 


PYTHON CODE : 








#importing the essential libraries 
from mpl_toolkits import mplot3d 
import scipy.special as spl 

import numpy as np 

import matplotlib.pyplot as plt 


import matplotlib.ticker as mtick 


#assigning the constant values for the cross section plot; 
w=124 

L=10 

g=8.8E-20 

x=0.2 


#defining our cross section function as z(x,y); 
def z_function(B, jy): 
return (np.sqrt(2*np.pi)*(Cy*L*g*B*(pow(x,2)))**2 * (1/16))*y*pow( 
x, 5)*np.exp(-(y**2)*(x¥**2))* 
spl.iv(0.5, (Cy¥**2)*(x**2))) 


#defining our variables- B=x and K=y; 
B = np.linspace(0, 1950, 100) 
y = np.linspace(0, 125, 100) 


#mapping the function z(x,y) to meshgrid Z(X,Y) 
X, Y = np.meshgrid(B, y) 
Z = z_function(X, Y) 


#plotting wireframe plot Z(X,Y) 

fig = plt.figure(figsize=(20,20)) 

ax = plt.axes(projection="3d") 

ax.plot_wireframe(X, Y, Z, cmap='winter') 

ax.view_init (32.5,-75) 

#labeling the axis for wireframe 

ax.set_xlabel('Magnetic field (B)'+ r'$\rightarrow$') 
ax.set_ylabel('Photon momentum (K)'+ r'$\rightarrow$') 
ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 
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A.1. 3D plots 


#fixing the figure size and zaxis format 
fig = plt.figure(figsize=(20,20)) 
ax = plt.axes(projection='3d') 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#plotting Z(X,Y) with color bar 

surf=ax.plot_surface(X, Y, Z, rstride=1, cstride=1, 
cmap='rainbow', edgecolor='none') 

fig.colorbar (surf , shrink=0.5,aspect=7) 

ax.view_init (32.5,-75) 

#labeling the axis for surface plot 

ax.set_xlabel('Magnetic field (B)'+ r'$\rightarrow$') 

ax.set_ylabel('Photon momentum (K)'+ r'$\rightarrow$') 

ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


#fixing the figure size and zaxis format 
fig = plt.figure(figsize=(20,20)) 
ax = plt.axes(projection="3d") 


plt.gca().zaxis.set_major_formatter (mtick.FormatStrFormatter('%.0e')) 


#plotting Z(X,Y) with color bar 

qcs=plt.contourf(X, Y, Z,20, cmap='rainbow') 
plt.colorbar(qcs,shrink=0.5,aspect=7) 

ax.view_init (32.5,-75) 

#labeling the axis for countour plot 

ax.set_xlabel('Magnetic field (B)'+ r'$\rightarrow$') 
ax.set_ylabel('Photon momentum (K)'+ r'$\rightarrow$') 
ax.set_zlabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 


#display our 3d plots 
plt.show() 
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.1. 3D plots 
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Figure 7.3(b) 


A.1. 3D plots 
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Figure 7.3(c) : contour. 
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A.2. 2D plots 


A.2 2D plots 


A.2.1 oi vs K plot 


PYTHON CODE : 








#importing the essential libraries 
from mpl_toolkits import mplot3d 
import scipy.special as spl 

import numpy as np 

import matplotlib.pyplot as plt 


import matplotlib.ticker as mtick 


#assigning the constant values for the cross section plot; 
B = 1950 

w=124 

L=10 

g=8.8E-20 

x = 0.2 


#defining our cross section function as z(x,y); 
def z_function(y): 
return (np.sqrt(2*np.pi)*(Cy*L*g*B*(pow(x,2)))**2 * (1/16))*y*pow( 
x, 5)*np.exp(-(Cy**2) *(x**2))* 
spl.iv(0.5, (y**2) *(x**2))) 


#defining our variables- R=x and K=y; 
y = np.linspace(0, 125, 50) 


z=z_function(y) 


#plotting wireframe plot Z(X,Y) 
fig = plt.figure(figsize=(16,10)) 
plt.plot(y,z,'r-o') 

plt.grid() 


#setting the axis format and range 
plt.gca().yaxis.set_major_formatter (mtick.FormatStrFormatter('%.1e')) 


plt.yticks(np.arange(0, max(z), 0.5E-33)) 


# Setting axis labels for the plot 
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A.2. 2D plots 


plt.ylabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$') 
plt.xlabel('Photon momentum (K)'+ r'$\rightarrow$') 


#display our 2D plots 
plt.show() 





OUTPUT : 
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Figure 7.4 : oto vs K plot . 
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A.2. 2D plots 


A.2.2 oi vs R plot 


PYTHON CODE : 








#importing the essential libraries 
from mpl_toolkits import mplot3d 
import scipy.special as spl 

import numpy as np 

import matplotlib.pyplot as plt 


import matplotlib.ticker as mtick 


#assigning the constant values for the cross section plot; 
B = 1950 

w=124 

L=10 

g=8.8E-20 

y=124 


#defining our cross section function as z(x,y); 
def z_function(x): 
return (np.sqrt(2*np.pi)*(Cy*L*g*B*(pow(x,2)))**2 * (1/16))*y*pow( 
x, 5)*np.exp(-(y**2) *(x**2))* 
spl.iv(0.5, (y**2) *(x**2))) 


#defining our variables- R=x and K=y; 
x = np.linspace(0O, 0.2, 50) 


z=z_function(x) 


#plotting wireframe plot Z(X,Y) 
fig = plt.figure(figsize=(16,10)) 
plt.plot(x,z,'r=o') 

plt.grid() 


#setting the axis format and range 
plt.gca().yaxis.set_major_formatter (mtick.FormatStrFormatter('%.1e')) 
plt.yticks(np.arange(0, max(z), 0.5E-33)) 

# Setting x axis label for the plot 

plt.ylabel('Cross section ($\u03C3_{tot}$)'+ r'$\rightarrow$ ') 


plt.xlabel('Solenoid Radius(R)'+ r'$\rightarrow$') 


#display our 2D plots 
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A.2. 2D plots 


plt.show() 
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Figure 7.5 : Oo: vs R plot . 
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APPENDIX B 


Acronyms and Abbreviations 


SM - Standard Model 
QCD - Quantum Chromodynamics 
P - Parity conjugation 
C - Charge Conjugation 
T - Time Reversal 
PQ - Peccei and Quinn 
WW - Weinbergand Wilczek 
BBN - Big bang nucleosynthesis 
WIMPs - Weakly Interacting Massive Particles 
CDM - Cold Dark Matter 
SSB - Spontaneous Symmetry Breaking 
KSVZ - Kim-Shifman-Vainshtein-Zakharov 
DFSZ - Dine-Fischler-Srednicki-Zhitnitshii 
eV - Electron Volt 
SNe la - Type Ia Supernova 
CMB - Cosmic microwave background 
CAST - CERN Axion Solar Telescope 
SUMICO - Tokyo Axion Helioscope Experiment 


ALP - Axion- Like Particles 


ie: 


PS-ALP - Psuedo Scalar ALP 
S-ALP - Scalar ALP 
RBF - Rochester-Brookhaven-Fermilab microwave cavity experiment 
UF - University of Florida 
ADMX - The Axion Dark Matter eXperiment 
HFET - Heterojunction Field Effect Transistor 
SQUID - Superconducting Quantum Interference Device 
DM - Dark Matter 
SN1987a - Type 2 Supernova 
LHC - Large Hadron Collider 
TPC - Time Projection Chamber 
RHESSI - Reuven Ramaty High Energy Solar Spectroscopic Imager 
HEMTs - High-Electron-Mobility Transistors 
BNL - Brookhaven National Laboratory 
RF - Radio Frequency 
IF - Intermediate Frequency 
FFT - Fast-Fourier transform 
LHe - Liquid Helium 
JPA - Josephson Parametric Amplifiers 
CAPSEr - Cosmic Axion Spin Precession Experiment 


NMR - Nuclear Magnetic Resonance 
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